Localized states at zigzag edges of bilayer graphene 
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We report the existence of zero energy surface states localized at zigzag edges of bilayer graphene. Working 
within the tight-binding approximation we derive the analytic solution for the wavefunctions of these peculiar 
surface states. It is shown that zero energy edge states in bilayer graphene can be divided into two families: 
(z) states living only on a single plane, equivalent to surface states in monolayer graphene; (if) states with 
finite amplitude over the two layers, with an enhanced penetration into the bulk. The bulk and surface (edge) 
electronic structure of bilayer graphene nanoribbons is also studied, both in the absence and in the presence of 
a bias voltage between planes. 

PACS numbers: 73.20.-r, 73.20.At, 73.21.Ac, 73.22.-f, 73.63.Bd, 81.05.Uw 



Introduction: The quest for new materials and material 
properties has recently led to graphene, the missing two- 
dimensional (2D) allotrope of carbon (]]]. Stability and bal- 
listic transport on the submicrometre scale, even at room- 
temperature, make graphene based electronics a promising 
possibility yj. Indeed, with Si-based technology approaching 
its limits, a truly 2D material with unconventional electronic 
properties is regarded with great expectations. 

Graphene is a zero-gap semiconductor, and this prevents 
standard logic applications where the presence of a finite gap 
is paramount. Band gaps can still be engineered by confin- 
ing graphene electrons in narrow ribbons Hi. However, 
the lateral confinement brings about the presence of edges, 
which in graphene can have profound consequences on elec- 
tronics. This is essentially due to the rather different behavior 
of the two possible (perfect) terminations in graphene: zigzag 
and armchair. While zigzag edges support localized states, 
armchair edges do not JfJ U 01 . These edge states occur at 
zero energy, the same as the Fermi level of undoped graphene, 
meaning that low energy properties may be substantially al- 
tered by their presence. The self-doping phenomenon |8; 
the edge magnetization with consequent gap opening 
among edge states driven effects. 

Bilayer graphene, as its single single layer counterpart, is 
also a zero gap semiconductor ifioll , but only in the absence 
of an external electric field: the electronic gap can be tuned 
externally Hill . Nevertheless, the question regarding the pres- 
ence of edge states in bilayer graphene is pertinent. Firstly, 
zigzag edges are among the possible terminations in bilayer 
graphene, and secondly, the presence of edges is unavoidable 
in tiny devices. 

In the present paper we show that zero energy edge states do 
exist at zigzag edges of bilayer graphene. An analytic solution 
for the wavefunction is given assuming a semi-infinite sys- 
tem and a first nearest-neighbor tight-binding model. The an- 
alytic solution we have found defines two types of edge states: 
monolayer edge states, with finite amplitude on a single plane; 
and bilayer edge states, with finite amplitude on both planes, 
and with an enhanced penetration into the bulk. A schematic 
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FIG. 1: Charge density representation for the two families of edge 
states in bilayer graphene at ka/2n = 0.35: monolayer (a), given by 
Eq. l !12b ; and bilayer (b), given by Eq. d 1 3 b - 



representation of the two families of edge states can be seen 
in Fig. Q] We also show that bilayer graphene nanoribbons 
with zigzag edges have four flat bands occurring at zero en- 
ergy, consequence of the two families of edge states localized 
on each edge. In the case of a biased ribbon, where the two 
planes are at different electrostatic potential and a band gap 
develops for the bulk electronic structure, the spectrum still 
shows two flat bands while the other two give rise to level 
crossing inside the gap. 

Surface states in semi-infinite bilayer graphene: The study 
of edge states in AB— stacked bilayer graphene given here 
is based on the ribbon geometry with zigzag edges shown in 
Fig. |2 We use labels 1 and 2 for the top and the bottom layers, 
respectively, and labels Ai and Bi for each of the two sublat- 
tices in layer i. Each four-atom unit cell (parallelograms in 
Fig. |2]) has integer indices m (longitudinal) and n (transverse) 
such that mai + n&2 is its position vector, where ai = a(l, 0) 
and &2 = — v3)/2 are the basis vectors and a w 2.46 A 
is the lattice constant. The simplest model one can write to de- 
scribe non-interacting electrons in AB-stacked bilayer is the 
first nearest-neighbor tight-binding model given by, 
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FIG. 2: Ribbon geometry with zigzag edges for bilayer graphene. 

where a^ mn (&i ;m n ) is the annihilation operator for the state 
in sublattice Ai (Bi), i = 1, 2, at position (m, n). The first 
term in Eq. (Q]i describes in-plane hopping while the sec- 
ond term parametrizes the inter-layer coupling (t±/t <C 1). 
Without loss of generality we assume that the ribbon has N 
unit cells in the transverse cross section (y direction) with 
n 6 {0, . . . , N — 1}, and we use periodic boundary condi- 
tions along the longitudinal direction (x direction). This last 
simplification enables the diagonalization of Hamiltonian (Q~|i 
with respect to the m index just by Fourier transform along 
the longitudinal direction, H = X^fc Hk, with given by, 

2 

/// = -tJ^XXiMK 1 + ^h^n + h;k,n-l] 
i—l n 

-t±y y j a\ ;kjn h-,k,n + h.c. . (2) 

n 

In order to search for zero-energy edge states we solve the 
Schrodinger equation, H k \fi, k) = E Ml k \fi, k), for E Mt k 
ii First we note that Hamiltonian Hk in Eq. (O effectively 
defines a ID problem in the transverse direction of the ribbon. 
It is then possible to write any eigenstate k) as a linear 
combination of the site amplitudes along the cross section, 
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Ia*, k) = ^2 2J [ a i( k > n ) K> k > n ) + n) \bi,k, n) ] , 

n i—l 

(3) 

where the four terms per n refer to the four atoms per unit 
cell, to which we associate the one-particle states |cj, fc, n) = 
c\. k |0), with Cj = a,, and i = 1, 2. In addition we require 
the following boundary conditions, 

ai(fc,iV)=a2(^AO=M^-l) = /32(^-l) = 0, (4) 

accounting for the finite width of the ribbon. Then, analo- 
gously to the single layer case, it can be straightforwardly 
shown that if Eq. © is a zero energy solution of the 
Schrodinger equation, the coefficients satisfy the following 
matrix equations: 
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where Dk = — 2cos(fca/2). As the 2x2 matrix in Eqs. © 
and © has the following property for any complex pk, 
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we conclude by induction that the general solution of Eqs. (0 
and © has the form: 
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for n > 1, where the matrix T„ is given by, 
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and T* is the matrix whose elements are the complex con- 
jugate of T n . One also requires the convergence condition 
|— 2 cos(/ca/2)| < 1, which guarantees that Eq. © is sat- 
isfied for semi-infinite systems. It is easy to see that the 
semi-infinite bilayer sheet has edge states for k in the region 
2n/3 < ka < 47r/3, which corresponds to 1/3 of the possi- 
ble fc's, as in the graphene sheet. The next question concerns 
the number of edge states. As any initialization vector is a 
linear combination of only two linearly independent vectors 
there are only two states per edge (per k). Moreover, Eqs (O 
and (O are edge states solutions on different sides of the rib- 
bon. When the semi-infinite system is considered only one of 
them survives. In particular, taking the limit N — > oo, and 
choosing the simplest linear independent initialization vectors 
[ai(k, 0), 0] and [0, a 2 (k, 0)], the two possible edge states are, 



and 



ai(fc,n) = ai(jfe,0)D^e _i ^ n 

a 2 (k,n) = -ai(fc,0)n7J^- 1 ^e-^("- 1 ) 



ai(k, n) = 

a 2 (k, n) = a 2 (k, 0)D^e" l ^ n 



(11) 



(12) 



Although linearly independent, it is clear that the edge 
states (fTTl i and (fl2l i are not orthogonal, except for ka = tt. 
It is convenient to orthogonalize (fTTT i with respect to (fl2l i so 
that we finally obtain, 

J ai(fc,7i) = a 1 (k,0)D™e- l *T n 

k (13) 

which, together with Eq. (fl2] i. represent all possible orthonor- 
malized zero-energy edge states for bilayer graphene. The 
normalization constants in Eqs. (fl2l) and ( fT3l are given by, 
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FIG. 3: (Color online) (a) - Charge density for bilayer edge states at 
ka/2n = 0.36. (b) - The same as in (a) at ka/2n — 0.34. 



An example of the charge density associated with Eq. (13[ 
is shown in Fig. @ for t±/t = 0.2, where the |ai(fc,n)| 2 
dependence can also be seen as the solution given by Eq. ( TTZb 
for \ot2(k, n)\ 2 , apart from a normalization factor. 

The solution given by Eq. ( fT2b is exactly the same as that 
found for a single graphene layer |5j,|g], where the only sites 
with non-zero amplitude belong to the A sublattice of layer 2, 
the one disconnected from the other layer. Solution dot , on 
the other hand, is an edge state that can only be found in 
bilayer graphene. The sites of non-vanishing amplitude for 
this edge state occur at sublattice A of layer 2, and at sub- 
lattice A of layer 1, which is connected to the other layer 
through t± (see Fig. Q]). Had we increased the semi-infinite 
sheet from the other side of the ribbon, and two similar edge 
states would have appeared in the opposite edge with non- 
zero amplitudes at sites of the B sublattices. In regards to the 
penetration depth, A, it is easily seen from Eqs. ( fT2l and ( [Oi l 
that both solutions have the same value: A = — l/ln|Z)fc|. 
Nevertheless, the solution given by Eq. ( fT3l l has a linear de- 
pendence in n which enhances its penetration into the bulk. 
We expect these states to contribute more to self doping then 
the usual single layer edge states @], as the induced Hartree 
potential which limits the charge transfer between the bulk 
and the edge will be weaker. Note that the key to self dop- 
ing is the presence of both an electron-hole asymmetry and 
extended defects. Electron-hole asymmetry may be due to in- 
plane next nearest-neighbor hopping (NNN) t', while edges 
play the role of extended defects. The finite t' shifts the en- 
ergy of edge states, leading to charge transfer between clean 
regions and defects. The energy shift for the single layer is 
given by Ek ~ —t'(Dl — 1) to first order in t', apart from a 
global factor of — 3t' 11211 . This is exactly the energy shift we 
get (away from the Dirac points) for bilayer graphene with in- 
plane NNN hopping, if we neglect terms of the order t't± /t 
and higher. 

Nanoribbons of bilayer graphene ( unbiased): So far we 
studied localized states at the semi-infinite bilayer graphene. 
Experimentally, however, the relevant situation is a bilayer 
ribbon. The band structure of a bilayer ribbon with zigzag 
edges is shown in Fig. [4] (a) for N = 400, obtained by nu- 
merically solving Eq. ©. We can see the partly flat bands 
at E = for k in the range 2tt/3 < ka < 47r/3, corre- 
sponding to four edge states, two per edge. The zoom shown 
in Fig. S] (b) for ka f» 27r/3 clearly shows that there are 
four flat bands. Strictly speaking, the edge states given by 
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FIG. 4: (Color online) (a) - Energy spectrum for a graphene bilayer 
ribbon with zigzag edges for TV = 400. (b) - Zoom of panel (a). The 
inter-layer coupling was set to t±/t = 0.2. 



Eqs. ( fT2b and dT3l l [and the other two resulting from Eq. (0] 
are eigenstates of the semi-infinite system only. In the ribbon 
the overlapping of four edge states leads to a slight dispersion 
and non-degeneracy. However, as long as the ribbon width is 
sufficiently large, this effect is only important at ka ~ 27r/3 
and ka ~ 47r/3 where A is large enough for the overlapping 
to be appreciable \7). As Eq. ( fT3l ) has a deeper penetration 
into the bulk, its degeneracy is lifted first, as can be seen in 
Fig. [4] (b). We may then conclude that edge states do ex- 
ist in bilayer graphene ribbons. We expect band gaps to open 
due to magnetic instabilities induced by electron-electron in- 
teractions, similarly to graphene single layer Si. Actually, 
the edge states we have found live only on a single sublattice: 
A or B depending on the edge they are localized. This kind 
of localized states favor a ferromagnetic arrangement along 
the edge and antiferromagnetic between edges 11311 . consis- 
tent with what is found by first principles for stacked graphitic 
strips Il4ll . Also half-metallicity should occur in graphene 
bilayer nanoribbons as a consequence of edge states, analo- 
gously to the single layer 11511 . 

From the point of view of scanning tunneling microscopy 
(STM) we notice that bilayer edge states give rise to different 
intensities depending on the ribbon edge. As an example we 
consider the ribbon shown in Fig. [2] and assume that the STM 
signal is essentially proportional to the local density of states 
of the upper layer. At the top zigzag edge the STM signal is 
due to edge states of the bilayer type, the only ones with finite 
amplitude on the upper layer [Eq. (1131) 1. On the other hand, at 
the bottom zigzag edge both bilayer and monolayer families 
have finite amplitude on the upper layer, and a higher STM 
intensity is expected therefrom. 

Nanoribbons of bilayer graphene (biased): It has recently 
been shown that the electronic gap of a graphene bilayer can 
be effectively controlled externally by applying a gate bias 
ill ill . We now consider the case of biased bilayer nanorib- 
bon with zigzag edges, where the presence of edge states 
should play a role. The bias gives rise to an electrostatic 
potential difference, V, between the two layers. This is 
parametrized by adding to the Hamiltonian in Eq. (Q3 the term 



— J2m,n( n l-- m n ~ ni-.m.n), With //,.,„ „ = 0,1. 



b\. m n bi-m,n- Edge states are strongly affected by the bias. 
The semi-infinite biased system has only one edge state given 
by Eq. (IT2l . as the edge state having finite amplitudes at both 



layers [Eq. < fl~3T >1 is no longer an eigenstate. In Fig. [5] we 
show the band structure of a bilayer ribbon for different val- 
ues of the bias. Two partially flat bands for k in the range 
2ir/3 < ka < Att/3 are clearly seen at E = ±V/2. These are 
bands of edge states localized at opposite ribbon sides, with 
finite amplitudes on a single layer [Eq. ( fT2b and its counter- 
part for the other edge]. Also evident is the presence of two 
dispersive bands crossing the gap. Both the closeness of these 
dispersive bands to E as ±V/2 for ka as tt and their cross- 
ing at E = near the Dirac points can be understood using 
perturbation theory in V/ 1. As surface states living at oppo- 
site edges have an exponentially small overlapping, and those 
belonging to the same edge are orthogonal, we can treat the 
solution given by Eq. cTT~3T > and its counterpart for the other 
edge separately. Starting with Eq. (fT~3T >. the first order energy 
shift induced by the applied bias is Ek = V/2((n*} — (n^)), 
where (n\) and (n*) give the probability of finding the local- 
ized electron in layer 1 and 2, respectively. The value of these 
quantities is easily obtained from Eq. ( TT31 through a real space 
summation, 
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The band dispersion is thus given by. 
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where the minus sign stands for the band of states localized at 
the opposite edge. The result of Eq. (fT8l is shown in Fig. [5] 
as a dashed line which is hardly distinguishable from the nu- 
merical result. Note that for ka m tt we have Dk — > 0, so 
that Et w ±V/2. This means that for ka ~ tt the edge state 
given by Eq. (fT~3T > is essentially localized at layer 1, which is 
clearly seen from Eqs. ( TToT l and (T7\ as long as t±/t <C 1. For 
1 — D 2 = t±/t the energy shift [Eq. (fT8l l is zero, which leads 
to band crossing. For t± <C t we can expand around the Dirac 
points, k^a = 27r/3, 4tt/3. If k = fco + dk, the crossing takes 
place for Ska = ±t±/ '(tv3), each sign being assigned to dif- 
ferent Dirac points. Note that 5k is V independent, and its 
value compares fairly well with the numerical results shown 
in Fig. [5] Indeed, the approximate result given by Eq. (fT8l 
only fails at the Dirac points, where the edge states localiza- 
tion Length diverges and their overlap has to be considered. 
Increasing the bias makes first order perturbation theory to 
break down. We have found numerically that the gap opens 
for V > t. 

Conclusions: We have shown that zero energy edge states 
do exist at zigzag edges of bilayer graphene. We have derived 
an analytic solution for the wavefunction assuming a semi- 
infinite system and a first nearest-neighbor tight-biding model. 
This analytic solution defines two types of edge states: mono- 
layer edge states, with finite amplitude over a single plane; and 
bilayer edge states, with finite amplitude over the two planes, 
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FIG. 5: (Color online) Energy spectrum for a bilayer ribbon with 
zigzag edges for different values of the bias V: (a) V = tx/10, (b) 
V — t±/2, (c) V = t±. Inter-layer coupling tj_/t = 0.2 and ribbon 
width TV = 400. The dashed lines are the analytical result [Eq. J18H . 



and with an enhanced penetration into the bulk. Edge states 
are present even in bilayer graphene nanoribbons, where edge 
magnetization as well as half-metallicity are expected to show 
up in analogy with single layer graphene. We have also shown 
the robustness of bilayer graphene edge states to the presence 
of an electrostatic potential difference between planes. 
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